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We have measured the spin fluctuations in the YBa2Cu30s.5 (YBCO6.5, T c =59 K) superconductor 
at high-energy transfers above ~ 100 meV. Within experimental error, the momentum dependence is 
isotropic at high-energies, similar to that measured in the insulator for two dimensional spin waves, 
and the dispersion extrapolates back to the incommensurate wave vector at the elastic position. 
This result contrasts with previous expectations based on measurements around 50 meV which were 
suggestive of a softening of the spin-wave velocity with increased hole doping. Unlike the insulator, 
we observe a significant reduction in the intensity of the spin excitations for energy transfers above ~ 
100 meV similar to that observed above ~ 200 meV in the YBCO6.35 (T c =18 K) superconductor as 
the spin waves approach the zone boundary. We attribute this high energy scale with a second gap 
and find agreement with measurements of the pseudogap in the cuprates associated with electronic 
anomalies along the antinodal positions. In addition, we observe a sharp peak at around 400 meV 
whose energy softens with increased hole doping. We discuss possible origins of this excitation 
including a hydrogen related molecular excitation and a transition of electronic states between d 
levels. 



PACS numbers: 74.72.-h, 75.25,+z, 75.40.Gb 

I. INTRODUCTION 

The behavior of the magnetic fluctuations of layered 
cuprates is strongly coupled to their superconductivity. 
The parent compounds are antiferromagnetic Mott insu- 
latorsi^ The materials are then superconducting if the 
CuC-2 planes are doped with holes whose concentration 
exceeds the critical value p c =0.055 (Ref. 3j) when the 
long range antiferromagnetism is largely suppressed. In 
the hole overdoped and non superconducting region, the 
superconducting phase is replaced by a Fermi liquid and 
the magnetic excitations are largely absent 4£ Despite a 
significant effort, there is no complete understanding of 
the nature of the magnetic excitations, or of their role in 
the superconductivity or their hole doping dependence. 
The magnetic excitations can be studied by neutron in- 
clastic scattering techniques but owing to the large ex- 
change constants resulting in excitations with energies 
well over ~ 100 meV— the measurements are difficult with 
reactor sources. Measurements at these high energies can 
be made using pulsed spallation sources such as the ISIS 
facility at the Rutherford Appleton Laboratory. 

The magnetic excitations of the parent undoped, an- 
tiferromagnetic, and insulating compounds have been 
well studied in both YBCO6.15 and La2CuC , 4» 7 -~— Both 
sets of experiments found the magnetic excitations to 
be dominated by well denned spin waves emanating 
from the commensurate antiferromagnetic Bragg peak at 



Q=(0.5,0.5) with a large velocity of he ~ 600 meV A with 
the top of the dispersion band located at ~ 300-350 meV. 
If the hole concentration is only slightly above the critical 
concentration, the magnetic scattering consists of short- 
range magnetic order and at low energies slow excitations 
characteristic of ordered magnetic regions separated by 
charge rich regions as we observed forYBCCVas (T c =18 

y\ 32.35 

For hole concentrations within the superconducting 
phase of YBCO, the low energy commensurate spin fluc- 
tuations are replaced by incommensurate rods of scatter- 
ing which disperse towards a commensurate resonance 
peak at a finite energy transferj 10 ' 11 The low-energy in- 
commensurate scattering is highly anisotropic with peaks 
only displaced along one reciprocal lattice direction ! 12 : 13 
The value of the incommensurate displacement scales 
with hole doping and with superconducting transition 
temperaturei 14 : 15 Both the hole doped YBCO6+2; and 
La 2 _a;Sr a ;Cu04 systems give the same values for the in- 
commensurability for a given hole doping 

The resonance peak defines a crossover energy that 
separates low-energy incommensurate fluctuations from 
high-energy fluctuations which have more symmetric line 
shapes in momentum than the low-energy incommen- 
surate fluctuations! 17 i 18 The resonance peak in neutron 
scattering appears as a commensurate peak in momen- 
tum and has a well-defined lincshapc in energy as mea- 
sured by a constant-Q scan at the antiferromagnetic 
zone center^ While there are theoretical and exper- 
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imental studies relating the resonance energy to su- 
perconductivity (Ref. [20I - I22T ) . the resonance energy 
does not universally scale with T c or the gap value 
where it has been observed, most notably in YBCOe+a; 
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(Ref. 19), La 2 -xSr x CuCU, Tl 2 Ba 2 Cu0 6+ a (Ref. 
Bi 2 Sr 2 CaCu208 + 5 (Ref. 13), and HgBa 2 Cu0 4 (Ref. 
For example, for YBCO6.5 the energy of the resonance is 
~ 33 meV while the superconducting temperature is 59 K 
and for Lai.875Bao.i 2 5Cu04 the resonance energy is found 
to be ~ 60 meV despite a superconducting transition 
temperature of less than 6 Ki 26 ' 27 This is further demon- 
strated by our measurements on YBCO6.35 (T c =18 K) 
where no undcrdamped resonance peak was observed. 32 
There has also been some debate on whether the inten- 
sity of the resonance peak is sufficient to account for the 
electronic anomalies ! 28 ' 29 While there is significant ex- 
perimental work correlating the resonance energy with 
the superconducting gap, these studies have typically fo- 
cussed on a narrow range of hole dopings near optimal 
hole concentrations. — We therefore conclude, based on 
these examples, that the resonance energy cannot rep- 
resent a universal magnetic energy scale connected with 
superconductivity in the cuprates. 

The momentum and energy dependence of the mag- 
netic excitations at energies higher than the resonance 
has also been a topic of some debate with reports suggest- 
ing that the momentum dependence resembles squares 
fRefe. l26l andl30h and others reporting circular patterns 
(Ref. 1311 ) consistent with spin-wave type of excitations. 
The energy dependence is also not well understood with 
some groups suggesting near vertical dispersion (Ref. l3fit ) 
and others suggesting nearly linear dispersion of the spin 
excitations. This behavior of the high-energy magnetic 
fluctuations is discussed in this paper. 

We have previously investigated the spin fluctua- 
tions in samples of YBCO6.35 (T c =18 K) and Ortho-II 
YBC0 6 .5 (T c =59 K)M& The magnetic excitation spec- 
trum has been mapped in very underdoped yet super- 
conducting YBCO6.35 (T c =18 K) over the entire band of 
magnetic excitations.— The low-energy fluctuations re- 
sult from short-range magnetic order and the slow mag- 
netic fluctuations are characteristic of regions of ordered 
spins separated by charge rich regions! 35 ' 36 The magnetic 
fluctuations disperse in a similar manner to the antifer- 
romagnetic insulator with commensurate spin-waves ex- 
tending up to ~ 200 meV. At these high-energies, the 
magnetic excitations become heavily damped in energy, 
display substantial spectral weight loss, and have a lower 
energy than that observed in insulating cuprates and pre- 
dicted from linear spin-wave theory. We proposed that 
this was due to the interaction of the magnetic excitations 
with particle-hole pairs above the so-called pseudogap en- 
ergy, for it matched the energy scales observed in photo- 
emission and thermal conductivity experiments! 37 ' 40 " — 
The pseudogap energy was found from the optical con- 
ductivity to decrease with increasing hole doping but its 
origin is not understood ! 38 ' 39 In this paper we show that 
the suppression of the neutron scattering cross section 



seen with neutrons also decreases with increased hole con- 
centration. 

We also made earlier measurements of the low- 
temperature magnetic response of superconducting 
Ortho-II YBCO6.5 up to ~ 100 meV. At energies above 
the resonance peak of 33 meV, the anisotropic low- 
energy incommensurate magnetic response is replaced 
by isotropic rings of scattering similar to the spin-wave 
response observed in the heavily underdoped supercon- 
ducting YBCO6.35 compound, but with a much reduced 
velocity^ The integrated intensity above the resonance 
energy agreed with linear spin-wave theory for two- 
dimensional magnetic fluctuations up to about 100 meV 
where the error bars in the measurements became large. 

Motivated by the high-energy results in YBCO6.35, in 
this paper, we have extended our investigation of the 
magnetic fluctuations in Ortho-II YBCO6.5 to higher en- 
ergies. We will show that the momentum dependence 
is isotropic in momentum at energy transfers above 100 
meV and, in contrast to our previous analysis the veloc- 
ity is larger and the excitations may originate from the 
incommensurate points. Furthermore we shall show that 
the intensity of the excitations shows a marked decrease 
at energy transfers of ~ 125 meV, an energy scale which 
agrees with the pseudogap energy for this particular hole 
concentration. 

Based on these results, we propose that the magnetic 
excitation spectrum of all cuprates are described by two 
energy scales, one representing the effects of the super- 
conducting gap and the other the pseudogap. A detailed 
comparison with other techniques will be presented. 

In addition to these results, we also present data from 
very high-energies beyond the top of the one-magnon 
band in the cuprates and find a mode with an energy of ~ 
400 meV. Its energy scales with the hole doping concen- 
tration and the momentum dependence is consistent with 
a magnetic origin. Various possible explanations for this 
peak are discussed including a direct orbital transition. 

The paper is divided into three main sections. The 
next section describes the experimental arrangement, the 
methods and the samples used. The following section 
describes the experimental results including a discussion 
of the momentum dependence, the dispersion and inte- 
grated intensity of the magnetic excitations. This sec- 
tion and an Appendix also discusses the mode observed 
at high energies near 400 meV. The final section consists 
of the summary and conclusions. 



II. EXPERIMENT 

We performed two different sets of experiments us- 
ing the MAPS and MARI direct geometry time-of-flight 
spectrometers at the ISIS facility of the Rutherford- 
Appleton laboratory. The MAPS position-sensitive- 
detectors are located 6 m from the sample and are de- 
signed so that there is no direction that has significantly 
poorer Q resolution arising from the shape of the de- 
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tectors. The incident beam was obtained via a Fermi 
chopper spinning at a frequency of 400 Hz with incident 
energies fixed to Ei=400 and 600 meV. The choice of 
these energies was based upon our experimental work in 
YBCO6.35 (T c =18 K) where the top of the magnetic scat- 
tering lay in the range of 200-300 meV. The sample was 
mounted with the [001] axis parallel to the incident beam 
in a closed cycle refrigerator. As described in detail in 
Ref. [U the value of L, and hence the bilayer structure 
factor, varies with energy transfer. By varying the angle 
of the [001] axis of the sample with respect to the incident 
beam, we showed that the optic and acoustic scattering 
intensities are equal at energy transfers above ~ 50-60 
meV in both YBCO6.5 and YBCO6.35. Since we are in- 
terested in excitations only above this energy we have not 
separated the acoustic and the optical contributions as we 
have done before for the low energy magnetic scattering 
around and below the resonance energy. Our results are 
independent of the value of L and we shall discuss the 
data only in terms of H and K . Since the [001] axis was 
parallel to the incident neutron beam, symmetry enables 
us to gain count rate by adding the four detector banks 
present at low-angles on the MAPS spectrometer. 

The experimental arrangement is similar to our previ- 
ous reports on YBCO6.5 with E.;=150 meV. The goal of 
this experiment was to measure to higher energies using 
Ei=400 meV and Ei=600 meV so that we could study 
how the magnetic excitations evolve close to the top of 
the one-magnon band. 

The spectrometer was calibrated by using a standard 
vanadium sample of known mass. The cross section was 
assumed to be isotropic in spin while the definitions 
°f x"(5') a; ) an d the total cross-section are the same as 
we used previously and agree with those of most other 
groups. The definitions of the magnetic cross section 
used by us are compared with those of other groups in 
Refs. |3l] andHil- We note that the vanadium calibration 
has been found to be in good agreement using an internal 
phonon standard. 

Motivated by the preliminary high-energy results ob- 
tained at the MAPS spectrometer, we continued investi- 
gations of the high-energy fluctuations using the MARI 
direct geometry time-of-flight spectrometer. MARI is 
particularly useful for studying the momentum depen- 
dence of an excitation as the detector bank extends to 
angles up to 26 ~ 130° whereas the MAPS detector bank 
only extends to 60° in scattering angle. This difference is 
particularly helpful to distinguish the magnetic scatter- 
ing from the scattering by the molecular modes involv- 
ing hydrogen^ The MARI detectors are 30 cm long and 
2.5 cm wide and located 4 m from the sample position 
and the measurements were conducted with the sample 
in vacuum without the presence of a cryostat or closed 
cycle refrigerator with the entire spectrometer at room 
temperature. This eliminated the possibility of observ- 
ing molecular modes from water condensed on the walls 
of the cooling system. For all the measurements we used 
the relaxed Fermi chopper to provide a monochromatic 



incident energy of 750 meV while spinning the chopper 
at 600 Hz. The sample was mounted with the [001] axis 
parallel to the incident beam. 

Both the MAPS and MARI measurements were made 
on a sample of Ortho-II YBCO6.5. The sample consisted 
of six orthorhombic crystals with a total volume of ~ 6 
cm aligned on a multi-crystal mount with a combined 
rocking curve width of ~ 2°. The samples were aligned 
such that reciprocal lattice reflections of the form (HHL) 
lay in the horizontal scattering plane. Details of the crys- 
tal growth, the detwinning by stress along the a axis, and 
the oxygen ordering have been published earlier 4^ The 
sample is partially detwinned with the majority domain 
occupying 70 % of the total volume. The oxygen or- 
dering had a correlation length of ~ 100 A in both the 
a and b directions. The DC magnetization shows that 
the superconducting onset temperature was 59 K with a 
width of ~ 2.5 K. The low temperature lattice constants 
were a=6. 81 A, &=3.86 A, and c=11.67 A. These crystals 
were used for our previously published studies^ Using 
MARI, we also studied a multi-crystal array of four ~ 1 
cc crystals of YBCO6.33 (T c =8 K) also coaligned in the 
(HHL) scattering plane. The YBCO6.33 materials were 
orthorhombic with lattice constants a=3.844 A, 6=3.870 
A, and c=11.791 A. The samples were twinned and no 
oxygen ordering Bragg peaks were observed with neutron 
scattering techniques. Since x rays showed orthorhombic 
symmetry over much shorter distances, only short oxygen 
chains were present. 



III. RESULTS 

In this section, we describe the results from our ex- 
periments at ISIS for the samples of YBCO6.5 (T c =59 
K) and YBC0 6 .33 (T c =8 K). The first three subsections 
focus on the Ortho-II sample of YBCO6.5 and describe 
our extended measurements for energy transfers up to ~ 
200 meV. We will discuss the dispersion of the magnetic 
excitations, the lineshape in momentum, and the inten- 
sity. In the final section we present our observation of 
a mode at ~ 370 meV in YBC0 6 .5 and at ~ 430 meV 
in YBCO6.33 where it illustrates that the energy of this 
mode scales with the doping. 



A. Dispersion - High velocity spin waves 
originating from incommensurate positions 

We performed constant energy cuts through the MAPS 
data taken with Ei=400 meV and 600 meV to extract the 
momentum dependence of the magnetic fluctuations as a 
function of the energy transfer. Representative cuts are 
presented in Fig. Q] with cuts taken along the [1 0] di- 
rection for Ei=400 meV and [1 -1] direction for Ei=600 
meV. The Ei=400 meV constant energy cuts shown by 
the left-hand panel capture the overall qualitative mo- 
mentum dependence of the magnetic fluctuations. At an 
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FIG. 1. Constant energy cuts through the correlated magnetic 
scattering with Ei=400 meV (left hand panels) and 600 meV 
(right hand panels) obtained with the MAPS direct geometry 
spectrometer. The cut directions were chosen to be [1 0] for 
Ei=400 meV and [lT] for E l= 600 meV. The Ei=400 meV 
data is integrated along K by ± 0.1 r.l.u. and the Ei=600 
meV is integrated along [-K,K] by ± 0.05 r.l.u. 



energy transfer of fi,w=63 ± 13 meV, the correlated mag- 
netic scattering near the antiferromagnetic wave vector 
position is well defined and shows little evidence for the 
presence of two displaced peaks. At higher energy trans- 
fers of 88 ± 13 meV and 113 ± 13 meV, the correlated 
magnetic scattering broadens considerably and develops 
a clear flat-top structure indicative of the dispersion of 
the magnetic excitations. The broadening and splitting 
becomes even more obvious at the higher energy trans- 
fers probed using Ej=600 meV and in particular for Tilo ~ 
163 ± 13 meV two distinct peaks are observed. The re- 
sults in Fig. Q] clearly show evidence for the dispersion of 
the magnetic excitations at high-energies above the res- 
onance. This finding differs from some experiments on 
YBC06.6 where vertical rods of scattering were proposed 
to exist above the resonance peak, and were ascribed to a 
Fermi surface explanation for the magnetic excitations. 30 

The dispersion of the magnetic excitations was quan- 
tified by fitting to two Gaussians symmetrically dis- 
placed from the antiferromagnetic position at Q=(0.5, 
0.5). This simple model describes the data well over a 
very broad energy range. The background was assumed 
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FIG. 2. The peak positions along the [1 0] direction, with 
respect to the (H,K) = (0.5, 0.5) position, of the magnetic re- 
sponse as a function of energy transfer for Ei=400 meV and 
600 meV. Previously published results taken with Ei=150 
meV are plotted for comparison illustrating the consistency 
of the data and the low-energy response. The solid line is the 
spin-wave velocity measured in YBCO6.35 of 7ic=590 meV A. 
The error bar in the velocity is ± 30 meV A. The dashed 
line originating from the incommensurate elastic position is 
discussed in the text. 



to be a constant plus a linear term with increasing mo- 
mentum transfer. The values of the peak positions in 
momentum as a function of the energy transfer are illus- 
trated in Fig. [5] with the momentum transfer along the 
[1 0] direction. The solid line in Fig. [5] shows the mag- 
netic excitations measured in YBCO6.35 (T c =18 K) with 
a slope of ?ic=590 meV A. The errorbar in the present 
measurement of the velocity is ± 30 meV A. The exper- 
imental data taken with E,=400 meV and 600 meV are 
represented by filled diamonds and circles. The previous 
data taken with Ei=150 meV are plotted as open squares 
and are in good agreement with the current experimental 
data set where there is overlap at fiu) ~ 75 meV. 

The magnetic fluctuations at energies much below the 
resonance have incommensurate wave vectors displaced 
along (0.5 ±S, 0.5), as described previously. These in- 
commensurate fluctuations disperse inwards and meet at 
the resonance at an energy of hto—33 meV. Above the 
resonance energy, the spin fluctuations disperse outwards 
and were originally interpreted by us as showing a signif- 



icantly reduced exchange constant or spin-wave velocity 
when considering data only below ~ 100 meV. Our new 
data, Fig. [2j taken with a larger incident energy extends 
the dynamic range of our previous results up to Tiu ~ 
200 meV. At high-energies the momentum of the mag- 
netic peaks agrees well with the previous measurements 
on heavily underdoped and superconducting compounds 
and insulating YBCOg+x namely linear spin-waves ris- 
ing from the commensurate antiferromagnetic position 
at Q=(0.5 0.5). The effective exchange constant at high- 
energies is larger and similar to the insulating value. 

An alternate interpretation of the data in Fig. [2] is 
that the high-energy wave vectors extrapolate back to 
the incommensurate wave vector at the Tilj=0 elastic po- 
sition. This interpretation is shown by the dashed line 
for which the spin-wave velocity, derived from fitting all 
the data above Tiuj = 40 meV, is given by ?ic=800 ± 50 
meV A. This value is about 1.3 times larger than the 
value of the velocity of ~ 600 meV A derived from the 
insulator. This velocity would correspond to an exchange 
constant of ~ 150 meV. While this number seems quite 
large compared with the parent YBCOg.o and La2Cu04 
compounds, it is similar to exchange constants derived 
in Nd 2 Cu0 4 and SrCu0 2 4i^ This model of the spin- 
excitations provides a good description of all of the data 
above ~ 50 meV. 

The extrapolation back to the incommensurate posi- 
tion is done to relate to the stripe model where the 
antiferromagnetic regions are separated from hole-rich 
regions 4^ It is then expected that cones of scattering 
should exist at low-energies below the commensurate 33 
meV resonance. However it has been shown in Ref. [5^ 
that for weakly coupled stripes the intensity of the in- 
ner branch is much larger than the outer branch. Such 
a stripe model could then provide an explanation for 
the inward type of dispersion observed here and in other 
cuprate materials. 



B. Pattern of spin fluctuations in momentum - 
isotropic planar rings above 33 meV 

We now discuss the detailed momentum dependence of 
the magnetic excitations at a fixed energy transfer. For 
two-dimensional spin waves, we would expect rings of 
scattering that are approximately isotropic in the a* — b* 
plane for a fixed constant energy slice. The shape of the 
scattering in momentum would change only at very high 
energy transfers near the top of the band where the spin- 
waves approach the zone boundary. Fig. [2] illustrates 
that we are far from this limit and well below the top of 
the expected one-magnon band. — ^ 

In contrast to spin- wave theory that gives isotropic cir- 
cles of scattering it has been suggested that the scattering 
could form a square in the case of ladder type excitations 
or when scattering originates from Fermi surface nest- 
ing^ir— The interpretation of the scattering in terms of 
squares has been proposed by several groups most no- 



5 



YBCQ6.5,E=138 +/-13meV 




02 04 06 08 1 

(H,0) (r.l.u.) 

FIG. 3. (Color online) A smoothed constant energy slice 
through the Ei=400 meV data at Hlj—138 ± 13 meV in the 
superconducting phase at 10 K. The dashed lines indicate 
the direction of one-dimensional cuts discussed in the text. 
The increased in intensity at larger momentum transfers (near 
Q=(l,l)) is the result of multiple phonon scattering present 
at higher energies. 



tably for La2- :c Ba. E Cu04 and also YBC06.6- We shall 
now consider which of these models describes our new 
data. 

A constant energy slice through the magnetic scatter- 
ing hoj = 138 ± 13 meV is displayed in Fig. [3] The 
data shows a ring of intensity centered at the antiferro- 
magnetic position of Q=(0.5, 0.5) plus a significant back- 
ground which increases steadily with momentum transfer 
Q and is the result of multi-phonon scattering. While 
causing a sloping background, the magnetic scattering 
above background is symmetric to within experimental 
around the Q=(0.5, 0.5) position consistent with a con- 
stant ring of scattering. 

To illustrate the isotropy, we display in Fig. @] one- 
dimensional cuts through the ring of scattering of Fig. 
[5] in two directions at 45° apart. The upper panel is 
a cut through the [1 0] direction while the lower panel 
is a cut through the [1 1] direction plotted as a func- 
tion of absolute momentum transfer (A -1 ) relative to 
the antiferromagnetic position. The sloping background 
has been subtracted to allow a direct comparison. As 
illustrated in Fig. [TJ a constant plus a slope provides 
an excellent description of the background over the en- 
tire energy range of these experiments. The solid lines 
are the sum of two Gaussians symmetrically displaced 
from the antiferromagnetic point. It can be seen that the 
data is well described by a model where the scattering is 
from isotropic spin- waves. Moreover, the spin momen- 
tum pattern in YBCO6.5 is inconsistent with a pattern 
of squares. If it were so for YBCO6.5, the peaks that oc- 
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YBCO Ortho-ll, Ei=400 meV 
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FIG. 4. One-dimensional cuts through the correlated mag- 
netic scattering displayed in the constant energy slice in Fig. 
[3] The data has been corrected for the sloping background 
described in the text. The upper panel is a cut along the 
[1 0] direction and the lower panel is a cut along the [1 1] 
direction. Both plots are displayed in terms of absolute mo- 
mentum transfer (A -1 ) and integrated over ± 0.08 A -1 and 
± 0.12 A -1 for [1 0] and [1 1] cuts respectively. 



cur at go=0.28 A^ 1 along the [1 0] direction would result, 
for the square pattern, in a peak at 0.40 A -1 along the 
[1 1] direction. This is not supported by our data and we 
conclude that the momentum pattern at high-energies is 
circular as expected for an isotropic cone of spin waves 
intersecting a plane of constant energy. The same con- 
clusion was obtained earlier based on more limited data 
around huj=78 ± 8 meV for YBCOg.5 and also for the 
magnetic excitations of heavily underdoped YBC06.35- 

This conclusion differs from that of several high- 
energy studies on the monolayer and bilayer hole-doped 
cuprates. For much more heavily doped YBCO6+2, Refs. 
[30 and [57l have reported a 45° rotation of the incommen- 
surability above the resonance energy with respect to the 
low-energy incommensurate peaks for YBCO. The low- 
energy incommensurability below the resonance along 
the a-direction in both samples is similar and corresponds 
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FIG. 5. The momentum-integrated intensity (per two Cu 2+ 
spins) is plotted as a function of energy transfer from mea- 
surements using both reactor and spallation neutron sources. 
The reactor and the Ei=150 meV data was previously pub- 
lished.—^ The current data taken with E 4 =400 meV and 600 
meV is illustated by filled circles and diamonds respectively. 
The data is integrated over the in-plane wave vectors along 
the a* and b* directions. 



to 5=0.1 while the incommensurability for our sample 
is (5=0.06 r.l.u. Since the incommensurability increases 
with doping, this shows that our samples are more un- 
derdoped than those materials where strong anisotropy 
has been observed. It is therefore possible that with in- 
creasing doping, the nature of the high-energy magnetic 
scattering changes. This is predicted by both Fermi sur- 
face nesting arguments and also stripe models (Ref. [H 
and l59h which show high-energy magnetic scattering to 
be very sensitive to the parameters of the theories. 



C. Intensity - Suppression of spectral weight above 
125 meV 



The spectrum of the local susceptibility for YBCO 6. 5 
(x"(cj) = J d 2 qx"(q, w)) from the magnetic scattering 
integrated over the two-dimensional wave vector is dis- 
played in Fig. [5] It summarizes all data we have taken 
at reactor (open circles) and at spallation sources (for 
Ei=150 meV, 400 meV, and 600 meV). The data was 
converted to absolute units through the use of a vana- 
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dium standard, as discussed previously, and is the abso- 
lute scattering cross section per formula unit of YBC0 6 5 . 
This is the scattering cross section for both bilayers and 
therefore for 2 Cu 2+ spins. We have not integrated over 
the third direction in wave vector because the scattering 
is assumed to be two-dimensional and the intensity forms 
a rod along the c* direction while well defined in the a*- 
b* plane. This assumption is valid at energies above ~ 
40 meV as demonstrated in Ref. [U The high-energy 
(Ej=400 meV and 600 meV) values were obtained by 
numerically integrating constant energy cuts similar to 
those represented in Fig. [TJ 

Spin-wave theory for a two-dimensional antiferromag- 
net gives a constant momentum integrated intensity at 
least in the region where the energy of the excitations 
varies linearly with wavevector. The proportionality 
constant requires the knowledge of the renormalization 
constant Z which has been measured and theoretically 
found to be ~ 0.5 for the case of the Cu 2+ spins in 
the cuprates^S This has been experimentally confirmed 
using neutron scattering for the insulating, parent anti- 
ferromagnets YBCOg.is and L^CuO^^ The expected 
constant value of the momentum-integrated cross section 
(as shown in Fig. [5]) is x"=3.5 y? B /eV for two Cu 2+ spins 
in the unit cell as shown by the dotted line in Fig. [SJ 

The momentum-integrated intensity at low-energies, 
below the 33 meV resonance peak, has been discussed ex- 
tensively 31 ! 33 While the resonance peak is strongly sup- 
pressed on entering the normal phase above 59 K, we 
have previously observed that there is a weak remanent 
peak present at this energy in the normal state at 85 K 
(T c =59 K). Given the uncertainty in the origin of the 
resonance peak and is relation to superconductivity, we 
focus our attention on the magnetic spectrum above the 
resonance energy. Data from MAPS using an incident en- 
ergy of 150 meV showed that the momentum-integrated 
intensity is roughly constant between the resonance en- 
ergy and below ~ 100 meV, and with the magnitude in 
agreement with spin wave theory. However, at high en- 
ergy transfers above ~ 100 meV, we now find that there 
is a marked and systematic decrease of the integrated 
intensity as shown by the data with Ei=400 meV and 
600 meV in Fig. [5] This is in direct contrast with linear 
spin-wave theory and the antiferromagnetic Mott insula- 
tor where the momentum-integrated intensity has been 
shown to be at least constant and then to peak at the 
top of the one-magnon band. 

This systematic suppression of fast spin fluctuations 
mirrors and reinforces our findings for the heavily under- 
doped YBC0 6 .35 (T c =18 K). There, we found that the 
high-energy fluctuations were damped in energy, lost sig- 
nificant spectral weight in comparison to the insulator, 
and were heavily renormalized to lower energies at the 
zone boundary. 34 This result was interpreted in terms 
of the spin-waves interacting with particle-hole pairs 
at high-energies above the pseudo-gap energy extracted 
from various transport measurements M- For YBCO6.35 
(T c =18 K), suppression of spin fluctuations occurred at 
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FIG. 6. A summary of the energy value at which the 
momentum-integrated susceptibility (x"( u )) fsHa to half the 
value predicted by spin-wave theory (taken to be 1.75 /i%/eV 
x 1/2 = 0.88 (i%/eV per copper spin) for a variety of hole 
doped cuprate materials. The light grey symbols are taken 
taken from a variety of techniques and summarized in Ref. 



around Tiu ~ 200 meV, while for more heavily doped 
YBCO6.5 (T c =59 K) suppression occurs at a much lower 
energy of ~ 100 meV. Because the upper energy scale de- 
creases as the doping increases we propose that the upper 
energy where the intensity suppression occurs follows the 
pseudogap. 

The loss of intensity at large energy transfers when the 
magnetic excitations are displaced significantly from the 
magnetic zone center has been observed in several metal- 
lic ferromagnets (such as Ni, Fe, and MnSi) i 61 i 62 In these 
systems the spin excitations become broad when both 
the momentum and energy transfer match those of the 
particle-hole pairs. A similar idea was also theoretically 
postulated to occur in lightly hole doped cuprates^ 3 - In 
contrast to the observation in magnetic metals, we ob- 
serve a disappearance of the excitations and not a redis- 
tribution of spectral weight. This result is particularly 
surprising as it is in violation of the total moment sum 
rule which states that the total integrated spectral weight 
must be a constant. It is possible that the spectral weight 
is very broad and shifted in energy resulting in it not be- 
ing distinguishable from background scattering. 

To confirm this point, and to directly compare our data 
with other cuprates, we have carried out an extensive 
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literature search of the neutron scattering experiments 
on hole doped cuprates where absolute normalization has 
been carried out and the energy value at which the high- 
energy integrated intensity (as plotted in Fig. [5]) falls to 
1/2 of the value predicted from linear spin- wave theory. 
We have taken this half height value to be 1.75 (j/g/eV 
per two Cu 2+ spins (or 0.88 ^? B /eV per Cu 2+ spin). 

Where there are inconsistencies for a particular hole 
concentration, we have taken the most recently published 
data for a given hole doped concentration. Most of the 
data has been taken with the MAPS spectrometer at 
ISIS, and so the data is at least self-consistent so far 
as the calibration to absolute units. We have previously 
compared (Ref. l3ll ) the absolute normalization for data 
taken on MAPS (using the vanadium standard) to an in- 
ternal calibration based on an acoustic phonon measured 
using a thermal triple-axis spectrometer and have found 
agreement within ± 20 %. The only exception is the 
data taken on nearly optimally doped YBCOg.ss which 
was taken with a triple-axis crystal spectrometer and no 
absolute calibration was performed. Here we have taken 
the energy to be 75 ± 25 meV based on the constant 
energy scans presented in the publication^ 

The data are plotted as a function of hole doping in 
Fig- El for hole doped cuprates YBa^usOe+z (Refs. loll 
and[3J)jLa2- a: Sr a ;Cu04 (Refs. lo5l - l69h . La 2 _ :r Ba :I ;Cu04 
(Refs. M and H3), and Bi 2 Sr 2 CaCu 2 8+(5 (Ref. [zO). 
The value of the hole do ping p was obtained from the 
formula proposed in Ref. l7ll . The light grey points are 
pseudogap values taken from a variety of techniques on 
a series of cuprates outlined in Ref. [lH and the refer- 
ences therein. The dotted line shows the curve proposed 
in Ref. [z3 (E pg = £ o (0.27 - x)/0.22) to describe the 
doping dependence of the pseudogap energy value. The 
trend of the magnetic data we present in Fig. [S]is in good 
agreement with this pseudogap formula with a value of 
E = 165 ± 10 meV for the maximum pseudoga p v alue in 
comparison to 152 ± 8 meV obtained in Ref. [72| to de- 
scribe a series of transport, optical, and tunneling data. 
The agreement between neutron and electronic probes is 
quite good especially given our pragmatic definition of 
the pseudogap as the value being the energy where the 
neutron integrated intensity falls to 1/2 the value pre- 
dicted by linear spin-wave theory. Our values for the 
trend of the pseudo-gap energy are consistent with those 
of Ref. [72| and enable us to conclude that the decrease in 
the high energy intensity suppression tracks and results 
from the decay of spin fluctuations to particle-hole pairs 
above the pseudogap energy. It is interesting to note that 
while the lower energy resonance only scales with dop- 
ing and concentration over a limited range near optimal 
doping, the high-energy scale tracks universally across all 
hole concentrations and cuprate concentrations. It is also 
consistent with probes directly sensitive to the electronic 
properties. 

Our experimental studies of the spin fluctuations there- 
fore indicate the presence of two distinct energy scales in 
the cuprate superconductors. The first energy scale is as- 



sociated with the resonance peak and scales with super- 
conductivity over a narrow range of doping and is not a 
universal energy scale across different cuprates and hole 
concentrations, as discussed in the introduction. Our 
high-energy work points to the presence of an additional 
energy scale which tracks the pseudogap energy and is 
consistent amongst different cuprate materials and hole 
concentrations. There is considerable experimental evi- 
dence for such a scenario from ARPES (Ref. |73T). Raman 
(Ref. [74]), and tunnelling measurements fRef. I75landl76l) 
and possible physical mechanisms have been discussed in 
Ref. [77] and [7|| The fact that we observe anomalous 
spin dynamics near the magnetic zone boundary may 
imply that this second energy scale is associated with 
electronic effects near the anti nodal positions. ARPES 
measurements have found evidence for such an energy 
scale which does not scale with the superconducting T c 
and decreases with increasing doping (Ref. [79t ). consis- 
tent with the present experimental results from neutron 
scattering. The energy scale was observed as a broad 
peak (termed the hump) along the (ir,0) direction which 
remains above the superconducting transition tempera- 
ture^ Therefore, we find the upper spin energy scale, 
observed here with neutrons, tracks the doping charge 
hump energy gap measured using ARPES. 



D. High energy mode - Detection, experimental 
concerns, and possible origins 

So far we have been concerned with the properties of 
the magnetic excitations within the band of the lowest 
energy excitation (one-magnon band). In this section, 
we examine the possibility of a magnetic response well 
above the upper energies associated with the one-magnon 
band at fruj ~ 250-300 meV. While this is a relatively un- 
reported region, there are several noteworthy x-ray and 
optical studies which reflect the need to explore this en- 
ergy range with neutrons. Measurements of the absorp- 
tion coefficient using mid-infrared spectroscopy have re- 
ported a strong peak at an energy of about 425 meV for 
insulating La 2 Cu04 and at a slightly lower energy of ~ 
375 meV for Sr 2 Cu0 2 Cl 2 i£i Relatively recently a peak 
was observed in La 2 Cu04 using resonant inelastic x-ray 
scattering at an energy transfer near 500 meVj 82 ' 83 While 
several explanations were proposed, the origin of this ex- 
citation is unclear. Since neutron scattering has both a 
well defined cross section and selection rules, it may help 
in identifying the origin of any magnetic peaks at these 
high energies. 

A search for these high-energy excitations was initially 
made using the MAPS spectrometer on YBCO6.5 and ev- 
idence was found for a peak at an energy transfer near 
about 370 meV. Due to the limited detector range on 
MAPS, we continued these measurements on MARI on 
samples of YBC0 6 . 5 (T c =59 K) and YBC0 6 . 33 (T c =8 
K). The results from the MARI experiment with an in- 
cident energy of 750 meV are plotted in Fig. [7] for both 
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FIG. 7. Constant Q cuts taken on the MARI direct geometry 
spectrometer with an incident energy of 750 meV and with 
the relaxed chopper spinning at 600 Hz. Panels a)-c) are for 
YBCOe.s (T c =59 K) and panels d)-f) for YBCOe.33 (T c =8 
K). The solid line is a fit to a gaussian function on a sloping 
background. 
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FIG. 8. The intensity of the high-energy peak measured on 
MARI for YBCOe.33 (T c =8 K, open circles) and YBCOe.s 
(T c =59 K, open squares) as a function of momentum transfer. 
The form factor from Ref. [13 is shown as a solid line. For 
comparison, the intensity as a function momentum transfer 
of the excited states of ZrH2 are also plotted to illustrate that 
the intensity is inconsistent with a hydrogen related molecular 
excitation. 



YBC0 6 . 5 and YBC0 6 . 33 . Both experiments on MARI 
were performed at room temperature without using a 
closed cycle refrigerator. 

A peak is confirmed for YBC0 6 . 5 (T c =59 K) at 370 
meV and is visible at low momentum transfer. As the 
momentum transfer increases the intensity decreases. A 
similar result is found for YBCO6.33 (T c =8 K), but with 
the peak shifted to a higher energy transfer of 430 meV 
and their energies plotted in Fig. [HI The solid lines in 
Fig. [7] are fits to a Gaussian on a sloping background. 

One reason for concern is that the peak energy is com- 
parable to vibrational frequencies from water or metallic 
hydrides that may be present in the sample as impu- 
ritiesi££~— Unfortunately, YBCOe+x may absorb small 
amounts of water. Given the large neutron cross section 
for hydrogen, even minute amounts of hydrogen could 
provide a significant cross-section which could be erro- 
neously interpreted as magnetic. The problem is even 
more complicated because the molecular modes involving 
hydrogen do not scatter at high-energies beyond a scat- 
tering angle of 20=90° because at large energy transfers 
the impulse approximation ensures that each scattering 



site behaves as an independent source. This approxima- 
tion has been reviewed and studied in detail in Ref. H3- 
As a result the standard rule that vibrational scattering 
grows continuously with momentum transfer while mag- 
netic scattering does not, must be treated with caution. 

We address this point in Fig. [8] where the intensity as 
a function of momentum transfer was extracted from the 
cuts shown in Fig. [7] and the two data sets for YBCO6.33 
and YBCO6.5 were normalized to agree for the lowest 
momentum transfer studied. The solid line is the calcu- 
lated Cu 2+ form factor from Ref. |84| and this agrees well 
with the data. 

The dashed lines are the results of a calculation us- 
ing harmonic theory for the momentum dependence of 
a molecular excitation involving hydrogen measured for 
ZrHo. 88 We have chosen ZrH2 as the experimental data 
have been extremely well characterized and modelled and 
also lower energy excitations (~ 150 meV) have been 
found to be quite similar in nature to those measured 
in YHj;. The n = 2 and n — 3 excitations plotted in Fig. 
[5] correspond to energies of 294 meV and 441 meV and 
therefore represent examples of the scattering from the 
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FIG. 9. The energy of the high energy mode is plotted as 
a function of hole doping for both LSCO and YBCO sys- 
tems. The YBCO inelastic neutron scattering (INS) data is 
described here, the YBCO mid-infrared was taken from Ref. 
[85l . the LSCO resonant inelastic x-ray scattering (RIXS) is 
from Ref. [8^ . and the LSCO mid-infrared results are from 
Ref. M. 



molecular excitations in this energy ranged We can see 
that whereas the momentum dependence of the molecu- 
lar excitations is peaked at a finite momentum transfer, 
for momentum transfers below ~ 15 A -1 , the intensity 
is inconsistent with our experimental data. 

Based on these results, there is considerable evidence 
against the excitations we observe resulting from molec- 
ular modes involving hydrogen. Firstly, the energy of the 
excitation changes by ~ 50 meV with a very small change 
in the oxygen concentration. Secondly, the momentum 
dependence is inconsistent with the previous measure- 
ments on metal hydrides at similar energies. Thirdly, if 
the excitation did originate from a molecular excitation 
involving hydrogen, then we would expect lower energy 
excitations which have stronger structure factors than 
the excitation at ~ 400 meV. A search revealed no such 
peaks. We therefore suggest that these results are mag- 
netic and are associated with scattering from the Cu 2+ 
ions. 

Despite these facts, we do note that a recent study on 
MARJ of La2Cu04 did observe a hydrogen related mode 
at ~ 450 meV which did not change energy with dop- 
ing.— While conclusive proof for the existence of an ex- 
citation will only come from corroborating evidence from 
other techniques, it is interesting to speculate as to the 
possible magnetic and electronic origins of this excita- 
tion. The energy of the peak approximately agrees with 
a variety of other techniques. The measured value of 430 



meV for heavily underdoped YBCO6.33 agrees quite well 
with the optical studies on the parent La2Cu04 com- 
pound in Ref. HH The energy is also in approximate 
agreement with a variety of resonant inelastic x-ray stud- 
ies that report a peak at 500 meV energy transferi 82 i 83 
The substantial decrease in the energy with hole doping 
has also been reported in Ref. |93| where the energy of the 
high-energy excitation was observed to decrease from 500 
meV to below 300 meV as the doping varied from the par- 
ent undoped compound to optimal concentrations. We 
summarize the energy positions measured by resonant 
x-rays, mid-infrared, and our neutron measurements for 
both LSCO and YBCO in Fig. EJ For YBC0 6+2; , the 
hole doping has been extracted from the value of T c us- 
ing the formula described in Ref. [tTJ- For consistency, 
we have omitted the data where only a Neel temperature 
exists except for the trivial case of p = 0. The trend with 
doping in Fig. [9] suggests that the peaks measured with 
neutrons and optical techniques have a common origin 
and also that there is a consistent softening with hole 
doping. We note that mid infrared experiments observed 
the peak to be very broad in energy while we observe a 
relatively sharp excitation in energy. 85 

Possible magnetic mechanisms for this peak include a 
crystal field (dd) type transition that might be expected 
from an orbital transition, or from a bi-magnon excita- 
tion suggested in papers on resonant inelastic x-ray scat- 
tering. The bi-magnon excitation is not the traditional 
two-magnon process investigated using neutron scatter- 
ing. It was shown in Ref. |86| that a finite magnon- 
magnon interaction can result in a resonant excitation 
at high-energies. It is interesting to note that the peak is 
predicted to occur at ~ 2.7 J which is ~ 370 meV, close 
to the energy observed here for YBCO6.5. 

An orbital magnetic transition is also a possible ori- 
gin for the high energy peak. We show in the Appendix 
that the local electrostatic field splits copper states into 
5 orbital singlets coupled to the spin. The result sug- 
gests that with the crystal field parameters arising from 
the local oxygen structure, the first excited orbital state 
could have the observed energy. Clearly, further studies 
are required to determine the underlying origin of this ex- 
citation and its relevance to the physics of the cuprates. 



IV. CONCLUSIONS AND SUMMARY 

We have investigated the high-energy spin response in 
the underdoped Ortho-II YBCO6.5 superconductor. Our 
most important finding is the presence of a second high- 
energy scale in the spin fluctuations which is character- 
ized by a suppression of spectral weight. A compari- 
son with techniques which directly probe the electronic 
properties illustrates that this second energy scale is as- 
sociated with the pseudogap. Within the error of the 
measurements, this second energy scale seems to scale 
universally across all cuprate compositions and dopings. 
Based on these findings, we find direct evidence for the 
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pseudogap in the spin fluctuations and evidence for the 
presence of two distinct energy scales in the cuprate su- 
perconductors. 

We have also searched for excitations above the top of 
the one-magnon band measured to be 300-350 meV in the 
parent La2CuC>4 compound. We have observed a well de- 
fined excitation at 370 meV in YBC0 6 . 5 (T c =59 K) and 
a similar excitation at 430 meV in YBCO6.33 (T c =8 K). 
The intensity as a function of momentum transfer is con- 
sistent with a magnetic excitation and not in agreement 
with expectations of molecular excitations involving hy- 
drogen known to exist in a similar energy range. This was 
demonstrated through a comparison to ZrH2, a typical 
and well characterized metal hydride. Excitations of sim- 
ilar energy have been observed in the mid-infrared region 
and also using resonant inelastic x-rays. While possibil- 
ities include an orbital dd transition and bimagnon type 
fluctuations, more work is required to understand the un- 
derlying mechanism and true nature of these excitations. 
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V. APPENDIX: CALCULATION OF SPIN AND 
ORBITAL EXCITATIONS IN YBCO e+I 

We have calculated the energies of high-energy transi- 
tions of the copper spin and orbital moments with the full 
symmetry of the local electric field environment and with 
spin-orbit coupling. We find that the magnetic dipole dd 
transitions occur close to the energies of peaks observed 
in neutron scattering experiments. 

It is commonly assumed that the magnetic properties 
of the Cu 2+ ion in cuprate superconductors may be de- 
scribed by a spin-one-half doublet in which the orbital 
magnetic moment is quenched. This assumption may 
give a good approximation for low-energy spin waves, 
being transitions between the two members of the spin 
doublet, but it fails to describe magnetic properties in 
which the orbital moment is active. The 9 electrons (or 
one hole) in the 3d shell of Cu 2+ form a Hund's rule state 
with spin S = 1/2 and orbital angular momentum L = 2 
coupled by the spin-orbit interaction 

H so — XL ■ S (1) 

where A=-102.9 meV for free ions<21 

The five L = 2 states, if subject to a local electric field 
having the cubic symmetry of an octahedron of neigh- 
bors, would be split into a ground doublet and excited 
triplet. A splitting pattern of similar symmetry occurs if 




O CM ' ,5 ' 



FIG. 10. (Color online) Local structure about a central copper 
ion (Cu21) in YBC06+;c. In the buckled a — b plane the 
neighbors of the Cu2\ atom (red) are four nearest neighbors 
oxygens, 02 and 03 (green) displaced up from the plane by 
dz, one apical oxygen (04'5' green) and four oxygens (02, 
03) in the bilayer plane above displaced down by dz. The 
n'th of the nine oxygen neighbors is denoted by the trailing 
label l n\ 

the interaction is regarded as arising from the hybridiza- 
tion of the Cu 2+ 3d states with the 2p-states of the O 2- 
neighbors. In the lower symmetry that prevails in YBCO 
one of the octahedral vertex neighbours in the z direc- 
tion is missing and the apical oxygen has only a minor 
effect on the Cu 2+ orbital ground state. The remaining 
four oxygen neighbors lie close to but are buckled out 
of the xy copper plane. This gives a strongly tetrag- 
onal field that splits the lowest cubic doublet by 338 
meV into the two orbital singlet states the lowest be- 
ing as 2 — y 2 = (| 2)+ I — 2))/y/2 as well as splitting levels 
of the upper triplet (sec Table I, left hand column). 

We have placed the neighbors at the low-symmetry po- 
sitions of the known crystal structured This gives the 
spectrum of excited states without the common assump- 
tion that a symmetric octahedral field dominates. The 
relative splitting pattern is expected to be reasonable 
while the overall magnitude of the splitting is expected 
to be less reliable. We regard the calculations as the 
lowest-order estimate of the energies of the excited or- 
bital states. 

The crystalline field energies have been calculated for 
up to 18 near neighbors but focus is on the 4 nearest 
2 ~ ions which are found to dominate the spectral pat- 
tern. The atomic structure around a central Cu 2+ ion, 
Cu2\ at (0,0,0), in the lower bilayer is shown in Fig. [TO] 
The buckling of the 021 and the 031 ions in the Cu{2) 
planes is included through their z fractional displacement 
dz— 0.019 as measured by Jorgensend The lattice con- 
stants are those of YBCO6.5. 

The method is described by Hutchingsd The local 
electric field contains terms up to fourth order in the 
electron coordinates. They are sums over neighbor sites 
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FIG. 11. Magnetic peaks in the dynamic susceptibility versus 
energy E for crystal field and spin orbit coupling in vanish- 
ing exchange field. S XX (E) gives in-plane fluctuations and 

Sneutron 

(E) gives their sum. The longitudinal fluctuations 
Szz (E) perpendicular to the plane of the orbital moment cou- 
pled to the spin are strong. Neutron resolution is 50 meV (full 
width at half height). The energy of the lower peak at 388 
meV lies close to that of the peak observed in our neutron ex- 
periments where a peak at 370 meV is observed. The neutron 
resolution was set to 100 meV. 

j of a geometric factor determined by local ionic posi- 
tions and proportional to the tesseral harmonics, Z%f(j), 
where m=2,4 and n=0...m. These are multiplied by an 
even, p = c, or odd, p — s order Stevens operator equiv- 
alent in L expanded as a symmetrized function of the 
non-commuting angular momentum operators, OJ^f(L). 
The site symmetry emerges because particular terms in a 
given Z^(j) coefficient summed over j cancel as required 
by symmetry. The matrix elements of the Stevens opera- 
tor equivalent yield the 5 by 5 matrix whose eigenvalues 
and eigenfunctions give the Cu 2+ energies and wave func- 
tions. One example is the 4'th order potential 

a nn n 

j 

C\ {Z\ c (j)0\ c + Z\ s {])0\ s ) + 
C\ (Zl c (j)0l c + Zf{j)Of) + 
C\{ZfO\ c + Zf{j)Ol s ) + 
Ci(Zt(j)Ot + Zr{3)Of)} (2) 

which enters the fourth order contribution to the Hamil- 
tonian as H4 — —e 2 94(r A )a^ x V4 where the Stevens factor 



FIG. 12. The dd magnetic transition spectrum for neutron 
scattering arising from crystal field, spin-orbit and exchange 
field of 250 meV. The largest peak arises from fluctuations 
along the c-axis S ZZ (E). S X x(E) gives in-plane fluctuations 
and Sneutron(E) gives their sum. The neutron resolution is 
50 meV. 

f?4=-2/63, and (r 4 ) = 2.498 is the mean fourth power of 
the 3d radial wave function in Bohr radii, arj. The simpler 
second order term (not shown) was calculated with 64=- 
2/21, and (r 2 ) = 1.028. The net crystal field Hamiltonian 
for the YBCO6.5 structure taking the oxygen charge to 
be Qj—-2 is 

H c = -104.90^ + 3.04O 2c - 1.09O^ + 

0.07Of - 14.10f (3) 

The coupling of the spin to the five orbital singlets yields 
a spectrum of ten L — S states. Each orbital singlet 
becomes a spin doublet. The spin-orbit interaction (Eqn. 
[1} serves to shift the energies of the purely orbital states. 
For example, the lowest two orbital singlets, for which 
the orbital moment is entirely off-diagonal and equal to 
L z =2, are split by 338 meV before spin-orbit coupling 
increases the splitting to 388 meV as shown in Table I. 

The crystal field and spin-orbit coupling were diago- 
nalized simultaneously. In Table 1, we show that the 
spectrum of Cu 2+ energies does not depend much on the 
number of neighbors. Henceforth, we present results only 
for the nn=4 O 2- in-plane buckled neighbors. 

The ground state wave function is then that of a pair 
of time-reversed states, which for spin up is 

|0,*) = 
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248 


247 


388 


390 


387 
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633 


601 


1168 


1175 


777 


1168 


1225 


800 


1188 


1314 


815 


1188 


1418 


1023 


1471 


1592 


1133 


1471 


1682 


1159 



YBCOs.s nn=9 2 ~ 
4 inplane, 1 (apical), 
4 (bilayer) 



H m/Z =250 meV 



TABLE I. shows the weak dependence of the lowest energy 
levels on the number of neighboring charges (nn). For 9 oxy- 
gen ions the results are similar to those with only the 4 nearest 
ions. For the exchange-free field the lowest purely orbital level 
at 338 meV is moved by the spin-orbit interaction to 388 meV. 
Here a=3.83, 6=3.86, c=11.78 A and the lower Cu02 plane 
lies at z =0.360. An exchange field in the ordered antifer- 
romagnetic would split the spin-generate ground degenerate 
state by 248 meV. 
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1471 











1471 


-0.02 


-0.47 






TABLE II. The table lists energies and gound-state transition 
amplitudes for crystal field and spin-orbit states in a vanishing 
exchange field. 



periment, in part because dp — 0.1 range of p is small 
compared with the nominal 2~ charge. We surmise that 
the peak energy declines with doping because of increased 
screening of the electric field by the growing carrier den- 
sity of the electron gas. To produce the very low peak 
observed with infra-red in Fig. [9]at ~ 200 meV would re- 
quire an unreasonable "effective" oxygen charge less than 
half nominal or a weakened hybridization. 

When magnetic order is present, the strength of the 
interatomic exchange is 



0.87 | 2, 1/2) - 0.01 | 0, 1/2) + 0.50 | -2, 1/2) 
+0.04|-l,-l/2). (4) 

This is very closely equivalent to | 0z) = 0.96 | x 2 — y 2 ) + 
0.26 | xy) which reveals the persistence of the strong 
x 2 -y 2 = (| 2)+ | -2))/V2 orbital contribution. The 
mixing of the only excited xy state occurs because the 
only strong spin-orbit effect on the ground state is via a 
single off-diagonal matrix element (xy \ L z \ x 2 — y 2 ) = 2 
coupling the orbital ground to excited states. 

The magnetic dipole spectrum is weighted by the 
squares of transition matrix elements between states. 
The neutron spectrum is shown in Fig. 1111 The peak at 
388 meV is a mixture of strong longitudinal zz symme- 
try fluctuations via |M Z | 2 and a spin lowering amplitude 
\M~\ 2 as listed in Table II. 

The energy of 388 meV calculated for the magnetic dd 
transition peak lies close to the observed energy of 370 
meV (Fig. [7]) in YBCO6.5. This inelastic peak arises from 
direct magnetic transitions as opposed to the bi-magnon 
mechanism adoped by some (Ref. llOOT ) to explain peaks 
seen in inelastic x ray scattering. Evidence for the higher 
energy peaks has been obtained with both midinfrared 
and resonant inelastic x ray scattering therefore provid- 
ing some validation of this calculationJ^r-iS^ 

To obtain the lower energies observed at larger dop- 
ing (Fig. [SJ we reduced the n ear-n eighbor oxygen charge 
to reflect hole doping (Ref. llOlh . This does lower the 
calculated peak energies but by less than observed in ex- 



H ex = jJ2Si-Sj (5) 

i,3 

where J=125 meV is greater than the spin-orbit interac- 
tion. For the ordered antiferromagnet a spin feels, in two 
dimensions, a static molecular held 

H mfz = US = 250meV. (6) 

The ground state wave function Eq. 2] for H m f z = is 
essentially unchanged at high-field. This is because an 
exchange held along z couples an up spin only to the up 
spin component of the distant excited off-diagonal orbital 
state. The result is that the ground state doublet splits 
symmetrically with negligible exchange mixing of excited 
states. 

A large moment, 2.01 /is, of which the orbital mo- 
ment is half, is produced when an exchange held H m j z 
of 250 meV acts along the c direction. The smallness of 
the increase over zero held (M z — 1.997) is because the 
z component of the orbital moment is off-diagonal be- 
tween the lowest two orbital states. The moment for a 
held in the plane is half as large. The uniform suscepti- 
bility is known to be larger along the c direction although 
NMR (Ref. Ogives a ratio of Xc/Xab of only 1.28 esti- 
mated from band-like anisotropy as the internal exchange 
and bilayer coupling of the real material serves largely to 
overcome the underlying single-ion easy axis. Our results 
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extend those of Abragam and Bleaney to relatively large 
spin orbit coupling, A/A = -0.30, where the xy state at 
Ao = 338 meV in the large tetragonal field of YBCO now 
lies below the 3z 2 — r 2 state of their perturbation theory. 
The susceptibility normal to the plane is again twice 
that within the planes with similar but reduced g- values 
(9z - 4 ). 

When the planar spin correlation length is long 
enough, a given spin may feel some effect of a root mean 
square local field allowing spin waves to propagate, as 
we have found not only for YBCO6.5 where low-energy 
incommensurate fluctuations extend over f - 20 A(Ref. 
[33h . but even more so for YBCO6.35 and YBCOs 33 where 
the elastic correlation lengths reach 42 A(Ref. l32land[35l) 
and 76 A(Ref. H3) respectively. This local effect we ig- 
nore for the highest energy transitions. 



E 


M" 


M+ 


M z 











2.01 


249 


1.89 








390 








-1.73 


633 


1.22 


-0.01 





1175 


1.43 


-0.18 





1225 


0.21 


1.03 





1314 











1418 








0.01 


1592 











1682 


-0.01 


-0.22 






TABLE III. The table lists energies and transition amplitudes 
for crystal field and spin-orbit states in an exchange field of 
250 meV. 

Since a mean exchange field of 250 meV has been in- 
cluded, the results of Fig. [12] simulate the single-ion 
spectrum of the ordered antiferromagnet. The calcula- 
tion with large exchange field along the z direction shows 
three peaks in the energy range of interest while the neu- 



tron spectrum of paramagnetic YBCO shows one, one 
that correponds to the calculated out-of-plane zz fluctu- 
ations. Existing measurements of antiferromagnetically 
ordered YBCO and LSCO have not to date detected such 
structure in the ~ 400 meV region, a difficult region for 
neutron spectroscopy. However, the calculated fluctua- 
tions have energies in doped YBCO comparable to our 
experiment. 

Table 111 shows the amplitude of the transition mag- 
netic momenta, M = L + 2S, whose square gives the 
magnetic peak intensity arising from transverse excita- 
tions of M + and M" and from out of plane fluctuations 
M z , which is longitudinal for exchange field in the z- 
direction. 

The "spin-wave" peak in Fig. [12] is the lowest energy 
peak and its intensity is proportional to the square of the 
transverse spin matrix element M x = (M + + M~)/2. 
Only M~ = 1.89 (Table III) connects the ground state to 
the first excited state making the spin-wave transition 
strength (M x ) 2 =0.89 while the longitudinal orbital tran- 
sition at 390 meV has a greater strength (M z ) 2 =2.99. 

Our calculations indicate that direct magnetic dipole 
transitions occur in the energy range accessible to neu- 
tron and x ray inelastic scattering and should not be ig- 
nored. This magnetic spectrum represents the insulating 
background upon which changes are made by hole dop- 
ing. Its spectrum is physically distinct from and of higher 
energy than the strong attenuation of magnetic spectral 
weight we have observed as the spin waves enter a pseudo- 
gap continuum above 250 meV. So me x ray papers have 
suggested (for example, Refill and llOdT ) that relatively 
low energies less than ~ 1 eV cannot arise from dd transi- 
tions. Our relatively simple estimates of spin and orbital 
magnetism suggest otherwise. The orbital moment is not 
quenched since it contributes half the ground state mo- 
ment, adds a substantial xy symmetry into the ground 
state, and gives a large intensity for inter-orbital tran- 
sitions. We have therefore shown that orbital moments 
play an important role in the energy spectrum when tran- 
sitions out of the ground state are made at a few hundred 
meV. 
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